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Self-organized criticality (SOC) observed in an abelian sandpile model (ASM) has been believed 
to be related to the absorbing phase transitions (APT) exhibited by a deterministic fixed energy 
sandpile (DFES). We critically investigate the fink between the SOC and the APT exhibited by 
the DFES. In contrast to the widespread belief, a phase transition in the DFES is shown not to 
be defined uniquely but to depend on initial conditions because of the non-ergodicity of the DFES. 
Furthermore, we show that a phase transition in the spreading dynamics of the DFES is intimately 
related to the percolation rather than the avalanche dynamics of the ASM. These results illustrate 
that the SOC exhibited by the ASM has nothing to do with phase transitions in the DFES. We 
discuss the implication of our result to stochastic models. 

PACS numbers: 05.70.Fh,05.07Ln, 64.70.qj,64.60.ah 
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The concept of self-organized criticality (SOC) coined 
by Bak, Tang, and Wiesenfeld (BTW) [ij has been 
invoked to describe the ubiquity of criticality in na- 
ture |2|j_l3j . Sandpile models with either determinis- 
tic [H 13] or stochastic Q toppling rule have played 
the role of Ising model in understanding SOC theoret- 
ically d, 0]- Although sandpile models were devised 
as prototypes with the spontaneous emergence of crit- 
icality without fine-tuning of any parameters, which is 
the characteristics of SOC [ij, the similarity of SOC to 
standard (nonequilibrium) phase transitions had been 
discussed immediately after the emergence of the con- 
cept The connection between SOC and standard 
phase transitions, especially absorbing phase transitions 
(APTs) III, has been materialized by defining a 
'fixed energy' ensemble of the sandpile model, or a fixed 
energy sandpile (FES) [IMi. The idea of the FES has 
elicited fruitful conclusions at least for stochastic sand- 
pile models; the FES version of the stochastic sandpile 
model is argued to form a universality class together 
with other APT models [13], which, in turn, has moti- 
vated to pave "paths to self-organized criticality" from 
conventional phase transitions 16, 17]. 

Despite the triumph of the idea of the FES, skepticism 
against the SOC-FES hnk has been unabated [el. [itI. [l8| . 
In comparison to the achievement in the analytic study 
of deterministic sandpile models, also known as abelian 
sandpile models ( ASMs) ^] , rigorous analysis in favor of 
the link between SOC and the FES still lacks. In this 



context, it is not surprising that the authors of Ref. [12| 
also worried about the possibility that the FES may ex- 
plore transient rather than recurrent (SOC) configura- 
tions. Thus, a first step to remove such skepticism would 
be to cement the link between the SOC in the ASM and 
its FES version or deterministic FES (DFES). The pur- 
pose of this work is to investigate if the above-mentioned 
program can be accomplished. 

However, this program does not seem fully successful. 



As an easy counterexample, it is clear that the avalanche 
dynamics of the one-dimensional BTW automata [l^ is 
completely different from dynamics at transition point of 
its FES version Although the ASM and the DFES 
in the one-dimensional chain do not share criticality, the 
SOC critical density is identical to the transition point 
of the DFES [H, Accordingly, it seems still valid 
to anticipate that the mean energy in an ASM hovers 
around the phase transition point Cc of the correspond- 
ing DFES by dissipation (when activity is too strong) 
and driving (when activity is absent) and that in the 
thermodynamic limit the mean energy C,s of the ASM at 
criticality is identical to C,c- Meanwhile, the ASM on a 
one-dimensional chain is rather pathological in the sense 
that the probability of observing any finite avalanche in 
the thermodynamic limit is zero [l9j . Hence, one might 
expect that except such pathological cases the program 
to associate the SOC with the DFES may turn out to be 
successful. 

Recently, however, even this expectation has been chal- 
lenged [l8|. By measuring in the thermodynamic 
limit which is preceded by the infinite time limit. Fey, 
Levine, and Wilson [l^ argued that is different from 
C^S for many different models. Because phase transitions 
in nonequilibrium systems are usually studied in the ther- 
modynamic limit followed by not preceded by the infinite 
time limit, the method employed in Ref. [18] might be 
considered inconsistent with the philosophy of the FES. 
Then can we restore the link between the SOC and the 
DFES if appropriate order of limits are taken? 

To answer the above question, let us begin with defin- 
ing the DFES. In the DFES on a graph with V ver- 
tices, each vertex starts with nonnegative integer energy 
Zi {i — 1, . . . ,V). A configuration change is governed 
hy aV X V integer matrix F which has properties such 



that < for i 9^ j, Y., 



F — 



(energy conserva- 



tion), and for any pair i,j of vertices there is a sequence, 
say a connecting sequence, of vertices with finite length 
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1 + 2 {ko = i,ki,..., ki+i = j) such that Fk^k^+i < 
for all < TO < ?. Except the energy conservation, these 
properties are same as in the ASM 6]. For convenience, 
z'[ = Fa — 1 is introduced and called the critical energy 
(at vertex i). We refer to a vertex i with energy higher 
than the critical energy as an unstable vertex and it can 
topple, resulting in the energy redistribution such that 
Zj — i> Zj — Fij for all j. It is clear that the total energy 
of a configuration C, E{C) = Ylii is conserved during 
toppling events. 

By a sequence of toppling events, a system may evolve 
to a configuration without an unstable vertex, which will 
be called an absorbing state. For convenience, we in- 
troduce a function of configurations A{C) which takes 
if the system starting from C will arrive at an absorbing 
state after a finite number of toppling events and 1 other- 
wise. Due to the abelian property of toppling events (3], 
this function does not depend on the way how toppling 
events are ordered and, accordingly, is well-defined for 
finite V . 

For the DFES, there are two bounds Zm and Zm such 
that 

Zm = min{E\\/c ^ithE{C)> E,A{C) = l}, (1) 
Z,„ = max{i;|Vc with E{C) < E,A{C) = 0}. (2) 

It is easy to get Zm — The existence of Zm is obvi- 

ous, but to find Zm for arbitrary F does not look simple. 
Fortunately, if Fij = 0, that is, if there is no greedy 
vertex Q, we can show that Z,„ = — 
J2i-^ii/'^- To this end, we notice that A{C) = 1 un- 
der the simultaneous parallel update (SPU) implies that 
the steady state is a limit cycle with a certain period 
T [H m. In one period, every vertex should topple 
exactly same number of times, say W, because there is 
no greedy vertex. Accordingly, the total energy flow be- 
tween vertices is 2WZm- Hence the total energy should 
not be smaller than Zm- 

To complete the proof, we will construct a configura- 
tion Cm with A{Cm) = 1 and E{Cm) — Zm as follows. 
Let P be a permutation of 1, . . . , V^. At first, the vertex 
P(l) = i is given energy Fa. Energy at other vertices 
are assigned by the following iteration; if the vertex P{k) 
has been assigned energy, the energy of vertex P{k + 1) 
becomes — Y^'j Fp(j,^i-) j, where the sum is over all ver- 
tices except . . . , P{k + 1). The finial configuration 
after all vertices are exhausted is Cm- It is obvious that 
E{Cm) — Zm- If an ASM is defined by the topphng 
matrix Aij — Fij + SijLi, where = 1 if the energy 
at vertex i of Cm is Fa and otherwise. Cm should be 
a recurrent state of this ASM. This is because by con- 
struction it should successfully go through the burning 
test [1]. Note that the burning test to find a recurrent 
state is applicable when there is no greedy vertex. If we 
add Li to all vertices to Cm and perform the toppling dy- 
namics according to A (with dissipation), the final state 



should be Cm after all vertices topple exactly once be- 
cause adding energy of amount Li to all vertices of a 
recurrent state of this ASM is an identity operation (3]. 
That is, by the DFES toppling rule. Cm should come 
back to Cm if the toppling order is taken exactly as the 
corresponding ASM, which implies A{Cm) = 1- Note 
that the minimum energy among recurrent states of this 
ASM is exactly Zm [Si]. Hence above consideration also 
proves that A{C) of a recurrent state C of a certain ASM 
constructed by making some vertices dissipative should 
be 1. 

Zm is in many cases larger than Zm, but the order 
relation is reversed for the DFES on a tree. By a tree 
is meant a system with \Fij\ = \Fji\ < 1 {i ^ j) and 
a unique connecting sequence between any two different 
vertices. In a tree, Zm is F — 2 and Zm isV—1. Hence 
for any configuration on a tree with E{C) > Zm = V 
AiO = 1. 

Since A{C) = for any C with E{C) < Zm and 
A{C) = 1 for any C with E{C) > Zm, one might expect 
that "the phase transition point" should lie between 
(m and C,m, where Cm (Cm) is the limiting value of Zm/V 
(Zm/V) in the thermodynamics limit V oo. Now we 
will investigate if this claim is valid. 

In the literature, the similarity between this phase 
transition and the usual APT has been emphasized fl3| . 
However, the DFES as an APT model has a certain prop- 
erty which is not possessed by usual APT models such 
as the contact process (CP) jl^l- In the CP, a finite 
system should fall into an absorbing (particle vacuum) 
state with probability 1 unless spontaneous death is pro- 
hibited. On the other hand, some finite systems in DFES 
cannot fall into an absorbing state. Since the phase tran- 
sition involves the thermodynamic limit as well as the 
infinite time limit, the order of these two limits is cru- 
cial. In the CP, if infinite time limit precedes the ther- 
modynamic limit, there is no nontrivial phase transition 
(everywhere inactive). However, if these two limits are 
taken in the opposite order, there is a nontrivial phase 
transition which is the main concern in the literature. In 
the DFES, however, nontrivial phase transition can be 
defined even if we take the infinite time limit first [lij . 
Our discussion about the APT exhibited by the DFES 
starts from this "unconventional" phase transition. 

To study an APT, it is also necessary to define the 
(mean) density of unstable vertices at stationarity, say 
p(C), where C is the initial configuration. If A{C) = 0, 
p{C) is unambiguously 0. However, p{C) depends on how 
the toppling events are ordered if A{C) = 1 [131 ■ For 
example, the devil's staircase [2l| is observable when the 
SPU is employed, but it should disappear if a random 
sequential update rule is employed because the steady 
state is not a limit cycle any more. In any case, it is 
trivially true that A{C) > p{C) for all C. To simphfy 
the discussion, we will exclusively assume the SPU. 

Let Po(C; C: ^) be the probability that an initial con- 
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figuration of the DFES with V vertices is C. Here, 
C is the mean energy per vertex satisfying V(^ = 
J2c E{C)Po{C; C, V). Let us introduce 

c 

h{(:;v) = Y,p{c)p^{c (3b) 

c 

and define "transition points" = 1 or 2) 

d')=irif{c lim 0,(C;1^)^O|. (4) 

Since (pi > 02, Cc^'' cannot be smaUer than c}^'^ 
41) 



By 



definition, C,rn < Cc' < Ca/- Although the equality c}c^ — 
Q seems plausible for any DFES, we could not find a 
general proof. We can prove this equality only for the 
DFES defined on a tree because the period of a limit 
cycle on a tree is either 1 or 2 [25| . which implies either 
A{C) = 2p{C) or A{C) = p{C) (that is, (/)2 < < 202)- 
Hence, c}c^ = ci^^ = 1 is satisfied on every tree regardless 
of the initial conditions. Note that the method used in 
Ref. [rS*] to find the transition point is equivalent to using 
01 with Poisson distributed initial condition. 

In fact, the definition (|4]) makes it impossible to define 
"the" transition point. Cc^' should depend on Pq- To ex- 
plain why, let us divide the whole configurations into two 
classes and 5^ such that if C e 5^ (C G 5'^ A{C) = 
(1). Due to the non-ergodicity of the DFES [21|, 5"^ and 
are disjoint and should exhaust all possible configu- 
rations. Let U{C:V) = T.c^s^MC;C.V)E{C)/V > 
(piZm/V. If one can choose Pq such that in the thermo- 
dynamic limit Ca /(C)^'(C - Co), where Cm < Co < Cm, 
9 is Heaviside step function, and /(C) which should be 
not larger than C is a certain function, this model shows 
a (possibly discontinuous) transition at Cc^^ = Co- Only 
in case Cm = Cm as in the DFES on a tree, the transition 
point dH) of the FES is unambiguously defined. Even 
worse, certain Po{C;(,V) can trigger reentrance behav- 
ior in the sense that the set oo) = 0} is not 
connected. In this context, the disagreement between Cc 
(actually Cc^^) and Cs observed in Ref. [3l is the generic 
feature of the DFES. 

To illustrate, we will investigate the DFES on a 
bracelet graph with periodic boundary conditions defined 
by Fij = ASij - 25|i_j|_i 
consideration given above, C 
Cc^'' should be located between 2 and 3. Now assume that 
the initial energy distribution at each vertex is drawn 
from 

p(2fc;C) = (l-Po(C))(| 



25|i_j| From the general 
- 2 and Ca/ — 3. Hence, 
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where fc is a non- negative integer and Po{C) is a certain 
function of C in the range < PoiO < 1- Due to the par- 
ity conservation at each vertex [l8| and its resemblance 
to the DFES on a chain, one can easily find that if C 
is smaller (larger) than 2 -l-po(C), '/'i(C;^) goes to zero 
(nonzero) in the thermodynamic limit. That is, phase 
transitions occur at points where 2 -|-po(C) — C changes 
sign. For example, let Po(C) = [4C + 1] (mod 2), where 
[x] means the largest integer not greater than x and the 
modulo 2 operation restricts the possible values to ei- 
ther or 1. Note that 2 + po{x) — x changes its sign 
at |,f,x- Hence, there is reentrance behavior in the 
phase diagram (absorbing, active, absorbing, then active 
phases as C increases). For other initial condition, one 
should use Po(C) = J2kLoP('^^ + ^' C) ^'^^ the criterion of 
phase transition. 

As argued in the beginning, the definition of transition 
points Eq. Q does not comply with the spirit of the usual 
APT. A proper definition consistent with the usual APT 



Cf)=irif{CI lim03(C,i)7^O}, (6) 
i3(C,0 - Jim y2n{Ct)Po{C-X,V), (7) 



where Ct is the configuration at time t from the initial 
configuration C, and n{Ct) is the density of unstable ver- 
tices averaged over time up to t [p{C) — limt^oo n(Ct)]- 



But note that C,^'^ cannot be larger than Cc^^ so it is 

(3) 

clear that Cc cannot be defined uniquely, either. 

Since Cm < Cs ^ Cm for a certain ASM constructed 
from the DFES with negligible dissipative vertices, it may 
be possible to find an initial condition which yields the 
criticality as in the ASM (fine-tuning of initial condi- 
tions). However, as we have seen, recurrent states are 
related to the active phase within our definition, so the 
criticality, if exists, has nothing to do with the recurrent 
states. If the ASM forms a universality class, the de- 
tailed structure of dissipative vertices should not yield 
a considerable difference. Hence, no DFES can exhibit 
the similar critical behavior as the ASM, if universality 
hypothesis is valid. 

The difference between the ASM and the DFES is even 
more striking when we consider a spreading dynamics, 
or a microscopic APT (mAPT) The mAPT stud- 

ies how an initially localized unstable vertices spreads 
through a background inactive region with infinite vol- 
ume (the density of unstable sites is zero for all time, 
hence the name "microscopic" ) . The order parameter in 
this case is the (survival) probability that the activity 
persists indefinitely. For illustration, our discussion will 
be restricted to the two-dimensional BTW model 1] . 

First, we construct the background inactive region in 
such a way that each site is assigned energy 0, 1, 2, or 3 
with the probability (l-p3)po, (l-p3)pi, (1-P3)P2, orp3. 



(2) 



p(2A; + l;C)=Po(C) 



C-1 



fe g-(C-l)/2 



k\ 



(5) 
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respectively (0 < Pi < 1 and po +Pi +P2 = !)• By vary- 
ing Pi's, one can study the mAPT with the background 
density ranging from to 3. After forming the back- 
ground density, we choose a site randomly (this is not 
a mathematically rigorous statement, but the rigor can 
be attained with ease) and perform the DFES dynam- 
ics. Obviously, no unstable site appears with probability 
1 — If p3 is larger than the site percolation thresh- 
old [27| and if the randomly chosen site is in the infinite 
cluster of sites with energy 3, the spreading dynamics 
never dies. Hence, the critical density Cmc is bounded by 



28| 



Cmc< (l-p*)(pi+2p2)+3p* 



(8) 



Hence by varying pi and p2 , one can have infinitely many 
transition points. In striking contrast to the previous 
discussion, the critical density can be smaller than — 
2. For example, if pi — P2 = 0, it is clear that = 
3p* ~ 1.778 24 < 2. Note that the numerical value of the 
site-percolation threshold is p* 



27| 



0.592 746 

Up to now, we have shown that phase transitions oc- 
curring in the DFES without greedy vertices have noth- 
ing to do with the SOC displayed by the corresponding 
ASM and that the program introduced in the beginning 
cannot be fulfilled. This naturally raises the following 
question. Can we make a similar conclusion for the re- 
lation between stochastic sandpile models and its FES 
version, or stochastic FES (SEES)? 

Clearly, one cannot apply the discussion in the above 
directly to stochastic sandpile models. Unlike the DFES, 
the mAPT of the SEES cannot be understood in the 
framework of the percolation. That is, even if there is a 
percolating cluster with low background energy density, 
it is very probable that this system will fall into an ab- 
sorbing state. Or even if there is no percolating cluster 
with sufficiently large background density, the stochas- 
ticity can allow the activity to jump to another (finite) 
cluster and can survive indefinitely with nonzero prob- 
ability. Besides, any finite system of the SEES, unless 
total energy is larger than Zm, should fall into an ab- 
sorbing state in the infinite time limit just as in the CP. 
In other words, a recurrent state of a finite system can 
evolve to an absorbing state of the corresponding SEES. 

However, the recent observation by Karmakar et 
al. ^28*1 suggests an indirect application of our discussion 
to the SEES. According to Karmakar et al. [2^, a certain 
deterministic sandpile model belongs to the same univer- 
sality class as the stochastic Manna model Since the 
critical density of a DFES cannot be uniquely defined as 
we have seen, a transition point of the FES version of 
the model in Ref. [28[ is not defined uniquely. If the uni- 
versality hypothesis in the sandpile models also implies 
the same universality class of the FES versions, the criti- 
cal behavior of the SEES would have nothing to do with 
the stochastic sandpile model. Of course, this claim be- 



comes meaningless, if the universality hypothesis is not 
applicable to the FES version or if the claim in Ref. 
is wrong. However, the numerical study in Ref. [2i 
rather convincing and the universality hypothesis seems 



valid [13|. In any case, we believe that it is still neces- 



sary to make the connection between the SEES and the 
stochastic sandpile models more firmly. 
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